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I. BASIC FACTS 

No free quarks are observed in Nature. The upper limit to 
the ratio of the abundance n q of quarks in ordinary matter to 
the abundance of nucleons n„ is — < 1(T 27 111 to be com- 

P Hp — 1 ' 

pared to the expectation in the Standard Cosmological Model 

^ « io- 12 H . 

The cross section for inclusive production of quarks in p-p 
collisions, a q is < 10~ 40 cm 2 (l| and should be of the order of 
the the total cross section i.e.w 10~ 25 cm 2 if quarks existed as 
free particles. 

An inhibition factor w 10~ 15 implies that the only natural 
explanation of what is observed is that n q and a q are exactly 
zero, protected by some symmetry. 

A possible deconfining transition is then a change of sym- 
metry, or an order-disorder transition, and can not be a cross- 
over. 

No experimental data exist on gluons . It is anyhow as- 
sumed that also gluons do not exist as free particles. 

Confinement is defined as absence of colored particles in 
asymptotic states. 

If QCD is the theory of strong interactions it should be able 
to account for confinement. 

In ref|3] the idea was first proposed that there could be a 
deconfining phase transition, at about the Hagedorn tempera- 
ture Th from the hadron phase to a plasma of quarks and glu- 
ons. The transition is being searched in heavy ion collisions at 
CERN SPS, at Brookhaven RHIC and will be studied at higher 
energies at LHC. No smoking-gun signal for the formation of 
quark gluon plasma is known. One of the main difficulties is 
to define operationally confined and deconfined. 

A deconfining transition has been observed in numerical 
simulations of QCD on the lattice, in the absence of dynamical 
quarks (quenched theory). 

In general the thermodynamics of a system is described by 
the partition function 



Tr[exp{-—)] 



(1) 



with H the hamiltonian . For a field system it can be proved 
that the partition function is equal to the euclidean Feynman 
integral with the time axis compactified from to y, and pe- 
riodic boundary conditions for boson fields, antiperiodic for 
fermions. 



d§exp[— I dl I d 3 xL(§(x,T))} 



(2) 



In lattice simulations this amounts to simulate the theory on a 
lattice NtNg with the space extension N s N t ,the time exten- 
sion. The temperature is then given by 



1 

aN, 



(3) 



with a the lattice spacing in physical units. Renormalization 
group at one loop gives the lattice spacing a in terms of the 
physical scale A^ of the lattice regularization scheme 



1 ^ 

a = -^exp\—) 



(4) 



with P 



2N 



the natural variable of lattice and bo 



"fjjtp'l-T'^ — the celebrated coefficient of the lowest 
order term of the p function, which is negative (asymptotic 
freedom). 

From eqs(3) and (4) we get 



A , P\ 
— exp{-— ) 

N t bo 
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N t exp{ M? ] 



(5) 



As a consequence of asymptotic freedom large T corresponds 
to small g ,or to order, small T to large values of g or to disor- 
der. In usual thermodynamic systems there is disorder at high 
T , order at small T. This unusual behavior naturally leads to 
the idea of duality 01 H|. 

Duality is a deep concept in statistical mechanics, string 
theory, field theory. It applies to systems in (d+1) dimensions 
with topologically non trivial spatial (d dimensional) configu- 
rations . The (1+1) dimensional Ising model is a prototype|5|. 
These systems have two equivalent complementary descrip- 
tions: 

1) A direct description in terms of local fields <P(x) , in 
which (<£>) are the order parameters . We shall generically 
denote by /j the topologically non trivial spatial non local ex- 
citations. This description is convenient in the weak coupling 
regime 

2) A dual description in which p are local fields, (jj) the 
(dis)order parameters and the fields <&(x) appear as non local 
excitations. The typical coupling of the dual go is w - and this 
description is convenient when go *C 1 or g ^ 1 , the strong 
coupling regime of the direct theory. 

The Ising model in (1+1) dimensions is defined on a square 
lattice in terms of a local field <t>(x) which assumes the val- 
ues ± 1 , with a ferromagnetic interaction between the nearest 
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neighbors, H = -J^$(/J).$(/f) and Z[P,4>] = exp(-ptf). 
The model is exactly solvable. It has a phase transition at 
some known p c : for p < p c the system is magnetized , i.e. 
(4>) ^ ; for p > P c (4>} = and the system is disordered. If 
one of the two dimensions no of the system is interpreted as 
time and the other one «ias space, the configurations with non 
trivial spatial topology are the kinks (antikinks) , which have 
spin <£> = — 1 for «i < x\ ( «i >x{) and +1 elsewhere at fixed 
«o- It can be shown that the operator /j which creates a kink 
also has values ±1 |5]|6], and that|5] 

Z[p,4>]=Z[p>] (6) 

with 

sinh ^=^m (7) 

The disordered phase in the direct description is ordered in 
the dual description. Duality maps the strong coupling disor- 
dered regime of the direct description into the ordered weak 
coupling regime of the dual description and viceversa. This 
suggests to look for the dual description of QCD in which the 
low temperature phase should be ordered , and to identify the 
topologically non trivial spatial configurations /j and the sym- 
metry of the dual description, i.e. the disorder parameters (/j). 

There are a number of candidate dual excitations , namely 
monopoles |7| , vortices|8], AdS5 strings 0. 

It is fair to say that the dual excitations are not yet known 
, but something can be said about the dual symmetry which 
should be responsible for confinement. 



II. THE STRATEGY 

The facts described above suggest a strategy for investigat- 
ing confinement in QCD : 

I) Check on the lattice if the deconfining transition is indeed 
order-disorder, i.e. check if there is a natural explanation in 
terms of symmetry of the huge inhibition to the existence of 
free quarks. Notice that the order of the deconfining phase 
transition can be studied without knowing the order parameter 
, by looking at scaling behavior of the specific heat , as we 
shall see below. 

II) Identify the symmetry responsible for confinement , and 
possibly the dual excitations. 

Of course one has first to detect a deconfining transition on 
the lattice. This task proves to be equally complicated as in 
nature. 

In the so called quenched approximation , or in pure gauge 
theory without dynamical quarks one can use the Polyakov 
line , which is the parallel transport in the time direction across 
the lattice. 

rr 

L(x) = TrPexp(i A°(x,t)dt) (8) 
Jo 

Due to periodic boundary conditions at finite T the path is 
closed and due to the trace L is gauge invariant. 



Consider now the large distance behavior of the space cor- 
relator of L 

D(x) = (L f (x)L(0)) (9) 

On the one hand D(x) is related to the static potential between 
a quark and an antiquark as 

V(x) = -TlnD(x) (10) 

On the other hand, by cluster property it behaves at large dis- 
tances as 

D(x)^ x ^\(L(0))\ 2 +Cexp(-^) (11) 

It follows that if (L) = 

y(x) w^ooOjc (12) 
which means confinement of quarks . If instead (L) ^ then 

V(x) ~ x ^oo constant (13) 

which means no confinement. On the lattice a transition is 
observed at T c w HQMev from a phase in which (L) = to a 
phase in which (L) ^ 0. The transition is first order for St/ (3) 
gauge theory, (L) is the order parameter, Z3 , the center of the 
group, is the symmetry. 

In the presence of dynamical quarks Z3 is not a symme- 
try, and therefore it can not be the symmetry responsible for 
confinement in nature. Nor can the string tension be an order 
parameter: even in the presence of confinement there is string 
breaking. When pulling apart a static q — q pair at some dis- 
tance the system prefers to convert the potential energy into 
quark pair creation. 

The situation for two flavor QCD is shown in Fig[l] , where 
the phase diagram is displayed. The two quarks are assumed 
for simplicity to have the same mass m 1 10]. 




FIG. 1: The phase diagram of Nf = 2QCD. The transition line is 
denned by the maxima of the specific heat and of the susceptibility 
of the chiral order parameter, m is the quark mass, /j the baryon 
chemical potential. 

At m w chiral symmetry, which is broken at low temper- 
ature, is restored at T c « HOMev . Above71- the chiral order 
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parameter vanishes. However at m ^ chiral symme- 
try is not a symmetry any more. The transition line in Fig[l] 
is defined by the maxima of various susceptibilities, the spe- 
cific heat, the susceptibility of the chiral order parameter, the 
susceptibility of the Polyakov line, which all coincide within 
errors. This coincidence means that when the line is crossed a 
rapid change occurs in the energy content, in the chiral order 
parameter, in the Polyakov line. Conventionally the region 
above the line is called deconfined phase, the region below 
it confined, but there is no way to test this statement oper- 
ationally. This is analogous to what happens for heavy ion 
experiments. The problem would be solved if there were an 
order parameter for confinement. 

An analysis can be made of the chiral transition (m = ) 
fTlll . assuming that the relevant degrees of freedom at the tran- 
sition are scalars and pseudoscalars. The order parameters are 
then the elements of the matrix: 



* = *y = (^ / (l+y s )>|fy> 



(14) 



By one loop renormalization group analysis of the most gen- 
eral lagrangean compatible with chiral symmetry, neglecting 
irrelevant terms , one looks for infrared stable fixed points, 
which are a necessary condition for a second order phase tran- 
sition. If no such fixed points exist the transition must be first 
order. If one of such points exists, the transition can be second 
order. For Nf < 3 no fixed point exists, so that the transition 
is first order, and this fits with lattice data| 12]. 

For Nf = 2 instead there are two possibilities: 1) The tran- 
sition is first order if the T|' mass my — > as T —>T C . Then 
the transition is first order also at m ;j u ^ and there is no tri- 
critical point at fj ^ 0. fA here is the chemical potential. (See 
Fig[l]) This possibility is consistent with the order disorder 
nature of the transition. 

2) The transition is second order. Then m^i ^ at T c , there 
is no transition at m ^ but only a crossover , and a tricritical 
point at /j ^ 0. The crossover is incompatible with an order 
disorder type of deconfining transition. It requires an unatural 
fine tuning in the theory as discussed in Section 1 . 

It is therefore of fundamental importance to determine the 
order of the transition by lattice simulations. 

The tool to do that is a renormalization group based tech- 
nique known as finite size scaling analysis. It consists in ana- 
lyzing the volume dependence of susceptibilities as the spatial 
size of the system L s is sent large, which is governed by the 
critical indexes. Explicitly for the specific heat Cv and for the 
susceptibility % of the order parameter one has 1 1 3 ] 1 10] 



C v -Co=L?<t>c(^,m^) 



and 



5C(w) - Xo = L] ® m) (xLj , mLf 
Here x = (1 — j-) ar, d 0C, 



(15) 



(16) 



y) and a, y 
indexes in the standard notation 



V and yt, are the relevant critical 
Their values for weak first 
order, second order 0(4) and 0(2) are listed in table I. 0(2) 
is included since the staggered formulation for the fermions 





yh 


V 


a 


Y 





0(4) 


2.487(3) 


0.748(14) 


-0.24(6) 


1.479(94) 


4.852(24) 


0(2) 


2.485(3) 


0.668(9) 


-0.005(7) 


1.317(38) 


4.826(12) 


MF 


9/4 


2/3 





1 


3 


V Order 


3 


1/3 


1 


1 


oo 



TABLE I: Critical exponents. 

which is used in the numerical simulations could break 0(4) 
to 0(2) at finite lattice spacing. 

The volume dependence of the susceptibilities measured on 
the lattice simulations is then compared to Eqs. (15) and (16) 
to determine the critical indexes, or the order and the univer- 
sality class of the transition. 

i 

There are two scaling variables in the problem, xL s v and 
mL y s h which makes the analysis difficult. Our strategy has been 
to keep one of the variables fixed in turn and to observe the 
dependence on the other one liToll . 

First we have varied m and L s by keeping mL] 1 ' fixed . To do 
that one has to assume a value for y/ t , which we have first put 
equal to that of 0(4) and 0(2) which happen to be equal [See 
table I] 1 10] and then equal to that of first order. The expected 
scaling is then 



and 



(C v -Co)/L? = F c (xZj) 



(X(w> - Xo)/LJ = F m (xZj ) 



(17) 



(18) 



The quantities on the left hand side of Eqs(17) , (18) should 
l 

only depend on xL s v . The scaling for 0(4) is shown in fig[2] 
. If there were scaling all the curves in the same figure should 
coincide within errors. The transition is not compatible with 
0(4). Similar results are obtained for 0(2). Preliminary data 
for the case in which mL] 1 ' is kept fixed assuming for yt, the 
value for first order 1 14] seem instead compatible with the cor- 
responding scaling equations Eqs(17),(18) 

Another possibility is to go large with L s keeping the vari- 
i 

able zLs fixed. Here large L s means that the extension of the 
lattice is much bigger than the inverse mass of the pion, say at 
least 10 times. The following scaling laws are then expected 



and 



(C v -Co)=m^Gc(zL]) 



X(w> ~ 5C° = mW G (w> ( xZ ^ ) 



(19) 



(20) 



These scaling laws can be checked both with 0(4)[0(2)] and 
with first order. The result is shown in Fig[3] and again fa- 
vors first order. Simulations in which approximations of the 
algorithm are eliminated and data are obtained with an exact 
algorithm 1 14] confirm these results. 

The problem is numerically very demanding, with CPU 
times of the order of Teraflops . year. We plan however to re- 
peat the analysis on larger lattices and with improved actions 
, since the issue "order-disorder" or crossover is fundamental. 
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FIG. 2: Check of Eqs.(17), (18) assuming 0(4) second order. No scaling observed. 




III. THE DUAL EXCITATIONS. THE DUAL SYMMETRY 

A popular candidate for dual excitations are vortices|8]. 
The number of vortices is a conserved quantity in (2+1) 
space time dimensions, not in (3+1) dimensions. More- 
over in the presence of quarks the symmetry Z3 is explic- 
itly broken. Gauge models with the same center have very 



different behavior at the deconfining transitions, and center- 
less systems like the gauge theory with gauge group G2 ex- 
hibit confinement! 18]. We shall then concentrate on mag- 
netically charged dual excitations (monopoles). A gauge in- 
variant field-strength tensor F^ v can be defined , the so called 
'tHooft tensor| 15], which is nothing but the field strength of 
the residual gauge symmetry after abelian projection. The in- 
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dex a, (a = 1 , . . , N — 1 ) refers to the (N- 1 ) different species of 
monopoles which exist in SU(N) gauge theory 1 7] 1161. The 
dual tensor is defined as usual 



F a = 

/J.V 



1 



pa 



and the magnetic current as 



(21) 



(22) 



The magnetic current is identically zero if Bianchi identities 
are satisfied . It can be different from zero in a compact for- 
mulation of the theory , in terms of parallel transports like the 
lattice formulation. Monopoles do indeed exist on the lattice. 
The magnetic current is in any case conserved due to the anti- 
symmetry of the 'tHooft tensor . 







(23) 



A topological symmetry not related by Noether's theorem to 
the form of the Lagrangean. This symmetry can either be ex- 
act (Wigner) , and then the magnetic charge is defined and su- 
perselected, or be Higgs broken. In that case the ground state 
has no definite magnetic charge and the system is a dual su- 
perconductor : magnetic charges condense like Cooper pairs 
(electric charges) do in an ordinary superconductor. This sym- 
metry is a good candidate symmetry for confinement . In the 
confined phase (T < T c ) the vacuum behaves as a dual super- 
conductor. The Chromoeletric field acting between a q — q 
pair is channeled by dual Meissner effect into Abrikosov flux 
tubes . Since their energy is proportional to the length of 
the tube, i.e. to the distance between the pair , this explains 
confinement. Above T c the Higgs phenomenon disappears 
, magnetic charge is superselected, and confinement disap- 
pears. Such a mechanism for confinement was proposed long 
ago in ref s 1 20 ] 1 2 1 ] . The idea has been directl y c hecked on 
the lattice. An operator ff{x) can be defined|[l3lO 13 • 
which creates a monopole of species a, (a = 1,..,N — 1) at x. 
In the deconfined phase one expects (p a ) —Q if there is su- 
perselection of magnetic charge, (p") ^ in the deconfined 
phase indicates Higgs breaking and dual superconductivity of 
the vacuum. 

The operators ^{x) are the lattice transcription of the con- 
tinuum operators 



H a (x,t) = exp[i / d 4 yE a ± (y,t)b ± (y-x)} 



(24) 



where b±(y — x) is the vector potential produced in y by a 
monopole sitting at x, £" is the color component of the chro- 
moelectric field along the direction of the residual gauge sym- 
metry which identifies the monopole 1 23] [ 16]. Since £" is the 
conjugate momentum to A" the operator fj a (x,t) simply adds 
to A a L the field of the monopole.lt is easy to show that 1 23 ] 



JdAexp[-$(S + AS)} 
JdAexp[-$S] 



It is then convenient to use, instead of (//') the quantity 



ap 



= (S)s-((S + AS)) ls+AS) 



(25) 



(26) 



Since Eq(25) implies that (p a ) 
gration of Eq(26) that 



1 at p = it follows by inte- 



(jj a ) = ex P [ P (pyp' 



(27) 



Lattice simulations show the following: 

1) For T < T c or p < p c , p tends to a finite limit as the spa- 
tial size of the lattice L s — > °°. In practice it becomes volume 
independent when the lattice size is larger than lfm. This im- 
plies by use of Eq(27) that (//'} ^ below the deconfining 
transition and hence that there is dual superconductivity 

2) For T > T c or p > p £ , i.e. above the deconfining transi- 
tion as L v — » oa 



-cL s - 



(28) 



with c > implying by use of Eq(27) that (//'} = in the ther- 
modynamic limit, or superselection of the magnetic charge. 
3) For T w T c or P w p c . the scaling law holds 



If L s m K ^> 1 the above formula becomes 

which implies the scaling law 

p" I 



(29) 



(30) 



(31) 



The scaling law Eq(31) can be used to determine the critical 
index v, and has the remarkable property of being m inde- 
pendent . It has been successfully checked for U(l) , SU(2), 
SU (3 ) pure gauge theories [17]|23] 1 25] giving the correct crit- 
ical index . For Nf — 2QCD it gives a very good scaling with 
the critical index corresponding to first order|27]. This is con- 
sistent with the results for the susceptibilities described in Sect 
2. 



IV. CONCLUSIONS 

We have argued that the only natural explanation of experi- 
mental data on confinement is that confinement is related to a 
symmetry, and therefore that the deconfining phase transition 
is an order disorder transition, and not a crossover. 

The primary goal is then to check this statement, which 
amounts to determine the order of the deconfining transi- 
tion. For that the study of specific heat is specially convenient. 
Preliminary results on Nf — 2 QCD around the chiral transi- 
tion are consistent with it. 

Only the existence of an order parameter allows to define 
operationally confined and deconfined . 

No definite statement can be made about the dual excita- 
tions of QCD. 

Something can be said on their symmetries : they carry 
magnetic charge. Vacuum is a dual superconductor in the 



???????????????????????? 



confined phase, and makes a transition to normal at decon- by simulations on larger lattices and with improved action, 
finement . 

At higher temperatures magnetic charge is superselected. This work was partially supported by M.I.U.R. Pro- 

We are working to refine and to assess the results reported getto:Teoria e fenomenologia delle interazioni fondamentali 
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